Interacting model of new agegraphic dark energy: observational constraints and age 

problem 
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Many dark energy models fail to pass the cosmic age test because of the old quasar APM 
08279+5255 at redshift z = 3.91, the ACDM model and holographic dark energy models being 
no exception. In this paper, we focus on the topic of age problem in the new agegraphic dark energy 
(NADE) model. We determine the age of the universe in the NADE model by fitting the observa- 
tional data, including type la supernovae (SNIa), baryon acoustic oscillations (BAO) and the cosmic 
microwave background (CMB). We find that the NADE model also faces the challenge of the age 
problem caused by the old quasar APM 08279+5255. In order to overcome such a difficulty, we 
consider the possible interaction between dark energy and dark matter. We show that this quasar 
can be successfully accommodated in the interacting new agegraphic dark energy (INADE) model 
at the 2cr level under the current observational constraints. 

PACS numbers: 95.36+x, 98.80.Es, 98.80.-k 



I. INTRODUCTION 

At its present stage, our universe is undergoing an ac- 
celerated expansion, which has been confirmed by many 
astronomical observations, such as, type la supernovae 
(SNIa) [1], large scale structure (LSS) [2] and the cosmic 
microwave background (CMB) [3], among others. All 
these observations indicate the existence of a "dark en- 
ergy" with negative pressure. The most important theo- 
retical candidate for dark energy is the cosmological con- 
stant, A, which fits the observations well, but is plagued 
with some severe theoretical difficulties, such as the so 
called "fine-tuning" and "cosmic coincidence" problems 
[4]. During the past decade, in order to unveil its nature, 
theorists have proposed many phenomenological mod- 
els of dark energy, e.g., quintessence [5], fc-essence [6], 
tachyon [7], phantom [8], quintom [9], braneworld [10], 
and Chaplygin gas [11]. Further, the possibility that dark 
energy might interact with dark matter, owing to their 
unknown nature, has also been seriously considered in 
many works to help solve the cosmic coincidence problem 
[12] and the cosmic doomsday problem [13]. For reviews 
of dark energy, see, e.g., [14]. 

In recent years, it has been found that many dark 
energy models get into trouble when tested by some 
old high redshift objects (OHROs). It is obvious that 
the universe cannot be younger than its constituents, 
so the age of some astronomical objects (at some red- 
shift), if measured accurately, can be used to test cos- 
mological models according to this simple age princi- 
ple. There have been some OHROs discovered, includ- 
ing, the 3.5 Gyr old galaxy LBDS 53W091 at redshift 
z = 1.55 [15] and the 4.0 Gyr old galaxy LBDS 53W069 
at redshift z = 1.43 [16]. In particular, the old quasar 
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APM 08279 + 5255 at redshift z = 3.91 is an important 
one, which has been used as a "cosmic clock" to constrain 
cosmological models. Its age is estimated to be 2.0 — 3.0 
Gyr [17]. These three OHROs at z = 1.43, 1.55 and 3.91 
have been used to test many dark energy models, includ- 
ing, the ACDM model [18], the general EoS dark energy 
model [19], the scalar-tensor quintessence model [20], the 
f(R) = yfR 2 - i? 2 , model [21], the DGP braneworld 
model [22], the power-law parameterized quintessence 
model [23], the Yang-Mills condensate model [24], the 
holographic dark energy model [25], and the agegraphic 
dark energy model [26]. These investigations show that 
the two OHROs at z = 1.43 and 1.55 can be easily ac- 
commodated in most dark energy models, whereas the 
OHRO at z = 3.91 cannot, even in the ACDM model [18] 
and the holographic dark energy model [25] . In this pa- 
per, we will investigate the cosmic age problem in the 
new agegraphic dark energy (NADE) model. We will 
show that the NADE model also faces the challenge of 
such an age problem. In order to escape from the cosmic 
age crisis, we consider the possible interaction between 
dark energy and dark matter. We will check whether the 
age problem can be solved successfully in the interacting 
new agegraphic dark energy (INADE) model. Of course, 
our discussions are based on the current observational 
constraints on the models. So, we will first place obser- 
vational constraints on the NADE and INADE models, 
and then discuss the cosmic age problem. 



II. THE NEW AGEGRAPHIC DARK ENERGY 
MODEL WITH INTERACTION 

In this section, we describe the INADE model in a flat 
universe. Many theorists believe that we cannot entirely 
understand the nature of dark energy before a complete 
theory of quantum gravity is established [27] . In the cir- 
cumstance that a full theory of quantum gravity is not yet 
available, it is more realistic to consider the possible cos- 



2 



mological consequences of some fundamental principles of 
quantum gravity. The holographic principle [28] is com- 
monly believed to be such a principle, so it is expected 
to play an essential role in investigating dark energy [29] . 
Along this line, a model of holographic dark energy has 
been proposed [30] (see also, e.g., [31] and [32]). The 
agegraphic dark energy model [33] is constructed in light 
of the Karolyhazy relation [34] and corresponding energy 
fluctuations of space-time. Actually, it has been proven 
that the agegraphic dark energy scenario is also a kind 
of holographic model of dark energy [33] . In such a holo- 
graphic model, the UV problem of dark energy is con- 
verted into an IR problem, since the dark energy density 
is inversely proportional to the square of the IR length 
scale, pde ~ L~ 2 . In the old version of the agegraphic 
dark energy model [33], the IR cutoff is chosen as the 
age of the universe, t (here it should be pointed out that 
the speed of light has already been taken to be 1, so time 
and length have the same dimension). However, there 
are some inner inconsistencies in this model; for details 
see [35]. Therefore, in this paper, we only discuss the 
new version of the agegraphic dark energy model [35] . 

In the NADE model, the IR cutoff is chosen to be the 
conformal age of the universe, 
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where n is a numerical parameter and m p is the reduced 
Planck mass. 

If we consider a spatially flat FRW universe containing 
agegraphic dark energy and dark matter, the correspond- 
ing Fricdmann equation is 

"^2 (P™ + Pi) ' ( 3 ) 



H 2 = 



3m; 



or equivalently, 



E{z) =m = (^\+f) 1/ \ (4 ) 



1 - Qg{z) 



where fl m o is the present fractional dark matter density 
and tt q = p q /(3m 2 H 2 ). Note that, for simplicity, in this 
paper we neglect the contributions from baryons and ra- 
diation. From (2), it is easy to find that 
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Obviously, fl m = p m /(3m 2 p H 2 ) = 1 - Cl q from (3). By 
using (1), (2), (3) and (5) and the energy conservation 
equation p m + 3Hp m = 0, we obtain the equation of 
motion for f2„: 



n' = fi q (i - n,) 3 — 




(6) 



where the prime denotes the derivative with respect to 
x = In a. Since 4- = — (1 + z) 4-, we get 



dQ q 

dz 



2(l + z) 



(7) 



From the energy conservation equation p q +3H(p q +p q ) = 
0, as well as (2) and (5), it is easy to find that the 
cquation-of-state (EoS) parameter of the NADE model 
is given by 
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The NADE model has been studied extensively; see, e.g., 
[36, 37]. In the following, we shall extend the NADE 
model by considering the interaction between dark en- 
ergy and dark matter. 

Assuming that dark energy and dark matter exchange 
energy through the interaction term Q, the continuity 
equations become 



p q + 3H (p q +p q ) = -Q 



and 



Pn 



3Hp m — Q. 



(9) 



(10) 



Owing to the lack of any knowledge of the micro-origin 
of the interaction, we simply follow other work on in- 
teracting dark energy and parameterize the interaction 
term generally as Q = 3H(ap q + /3p m ), where a and f3 
are dimcnsionlcss coupling constants. In order to reduce 
both complexity and the number of parameters, one often 
considers the following three cases: (i) a = b and /3 = 0, 
denoted as INADE1, (ii) a = and j3 = b, denoted as 
INADE2, and (iii) a = f3 = b, denoted as INADE3. Note 
that according to our convention b > means that dark 
energy decays to dark matter, while b < means that 
dark matter decays into dark energy. In cases (i) and 
(iii), a negative value of b would lead to the unphysical 
consequence that p m becomes negative in the distant fu- 
ture. For negative b values in case (ii), no such difficulty 
exists. In [38] it is argued from the thermodynamical 
view that the second law of thermodynamics strongly fa- 
vors that dark energy decays into dark matter. So, in 
general, b is taken to be positive. 

However, recently, it has been found that observations 
may favor the case of dark matter decaying into dark en- 
ergy [39, 40]. In particular, in [41], in a way independent 
of specific interacting forms, the authors fitted for the in- 
teraction term Q using observations. They found that Q 
is likely to cross the non- interacting line (Q = 0), i.e., the 
sign of the interaction Q is changed, around z = 0.5. This 
raises a remarkable challenge to the interacting models, 
since the general phenomenological forms of interaction, 
as shown in the above, do not provide the possibility of 
changing signs. As noted in [41], more general forms of 
interaction should be considered. For this reason, a new 
form of interaction, (iv) a = —j3 = b, denoted as IN- 
ADE4, was considered in [42]. Obviously, for this case, 
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in the early stage, since p m > p q , Q is negative. How- 
ever, Q may change from negative to positive when the 
expansion of the universe changes from decelerating to 
accelerating. The parameter b in this case is also as- 
sumed to be positive, since negative b would lead to a 
negative p m in the distant future. 

For clarity, we denote the aforementioned intcration 
term Q as 



Q 



3bHp q , 
3bHp m , 
3bH(p q + p m ), 
3bH(p q - p m ). 



(11) 



Note that in cases (i), (iii) and (iv) the parameter b is 
always assumed to be positive in the literature. However, 
in the present work, instead of making such an assump- 
tion, we permit this variable to be completely free and 
allow the observational data to tell us the true story, no 
matter if the ultimate fate of the universe is ridiculous or 
not. 

Differentiating (5) with respect to In a and using (1), 
we get 



n' = n a l -2 



H 



^0 



(12) 



Differentiating (3) with respect to time t and combining 
(1), (5), (9) and (10), we can easily hnd that 
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Therefore, we obtain the equation of motion for D, q , 
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Qi 



(14) 



or equivalcntly, 

dz 1 + z 
where 



3 - 2(1 + z) ya 



Qi = 



Q 



3m 2 H 3 



(15) 



(16) 



From (5) and (9), we get the EoS parameter of dark 
energy: 



1 + ^V^9-Q2, 

dna 



where 
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3H Pq 
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It is convenient to define the effective EoS parameters for 
dark energy and dark matter as 
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and 
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According to the definition of the effective EoS parame- 
ters, the continuity equations for dark energy and dark 
matter can be re-expressed in the form of energy conser- 
vation: 
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3H(l + w^)p q = 



and 
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Taking the aforementioned four cases of interaction, one 
can obtain 



-b 



i-n„ 



= < 



-b 
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Q = 3bHp q , 
Q = 'SbHp m , 
Q = 3bH(p q + p m ), 

Q = 3bH(p q - p m ). 



(23) 



Now, the Friedmann equation can be expressed as 

H(z) = HqE(z), (24) 

where 

1/2 



E(z) 



{i-n q0 )e 3 fo 
i-n q {z) 



(25) 



In the above equation, Q q (z) can be obtained by numer- 
ically solving (15) with the initial condition VL q (zi n i) = 
?i 2 (l + z ini )- 2 /A at z lm = 2000 [36]. While this initial 
condition is obtained from a NADE model without inter- 
action in the matter-dominated epoch, it is still suitable 
to an interacting model of NADE because the contri- 
bution of dark energy to the cosmological evolution is 
negligible in the matter-dominated epoch and therefore 
its impact can be ignored during that time; for details 
see [43]. 



III. OBSERVATIONAL CONSTRAINTS 

In this section, we place observational constraints on 
the NADE and INADE models. For the data, we will use 
the combination of SNIa, CMB and BAO. 

First, we consider the 557 Union2 SNIa compiled in 
[44]. The theoretical distance modulus is defined as 



Pth(zi) = 51og 10 D L (zi) + po, 



(26) 
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where fj,Q = 42.38 — 5 log 10 h with h the Hubble constant 
H in units of 100 km s _1 Mpc -1 , and the Hubble-free 
luminosity distance 



D L {z) = (1 + 2) 



dz 



£7(2; p)' 



(27) 



where E = H/Hq, and p denotes the model parame- 
ters. Correspondingly, the \ function for the 557 Union2 
SNIa data is given by 



1=1 



a 2 (zi) 



(28) 



where a is the la error in distance modulus for each 
supernova. The parameter [1q is a nuisance parameter 
but it is independent of the data points. Following [45], 
the minimization with respect to /io is easily performed 
by expanding \ 2 of (28) with respect to /io as 



xl(p) = A-2 f i B + f i 2 C, 



(29) 



where 
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and 
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It is evident that (29) has a minimum for /io = -B/C at 



C 



(30) 



Since xl, mm = xl, instead of minimizing xl we will 
minimize x 2 which is independent of the nuisance param- 
eter /io. Obviously, the best- fit value of h can be given 
by the corresponding /i = B/C of the best fit. 

Next, we consider the cosmological observational data 
from WMAP and SDSS. For the WMAP data, we use 
the CMB shift parameter R; for the SDSS data, we 
use the parameter A of the BAO measurement. It is 
widely believed that both R and A are nearly model- 
independent and contain essential information of the full 
WMAP CMB and SDSS BAO data [46]. 

The shift parameter R is given by [46, 47] 



R = Q 



1/2 



dz 

W) 



(31) 



where the redshift of recombination = 1091.3 which 
has been updated in the WMAP 7-year data [48]. The 



shift parameter R relates the angular diameter distance 
to the last scattering surface, the comoving size of the 
sound horizon at z* and the angular scale of the first 
acoustic peak in the CMB power spectrum of temper- 
ature fluctuations [46, 47]. The value of R has been 
updated to 1.725 ± 0.018 from the WMAP 7-year data 
[48]. On the other hand, the distance parameter A from 
the measurement of the BAO peak in the distribution of 
SDSS luminous red galaxies [49] is given by 



.4 
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dz 



£7(2) 



2/3 



(32) 



where zt, = 0.35. In [50], the value of A has been de- 
termined to be 0.469 (n s /0.98)-°- 35 ± 0.017. Here, the 
scalar spectral index n s is taken to be 0.963, which has 
been updated from the WMAP 7- year data [48]. So the 
total x 2 is given by 



X' 



~ 2 2 2 

Xa + Xcmb + XbAOi 



(33) 



where xl is given in (30), Xcmb = ( r ~ Robs) 2 / o 2 r and 
Xbao = (A — A b s ) 2 /a A . The best fit model parameters 
are determined by minimizing the total x 2 ■ The 68.3% 



confidence level is determined by A% 



X Xn 



< 1.0, 



2.3 and 3.53 for 



1, 2 and 3, respectively, where 



n p is the number of free model parameters. Similarly, 
the 95.4% confidence level is determined by A% 2 = 
X 2 ~ Xmin < 4 A 6 - 17 and 8 - 02 for n p = 1, 2 and 3, 
respectively. 

Now, let us discuss the observational constraints on the 
NADE model. The NADE model is a single-parameter 
model whose sole parameter is n. Solving (7) numerically 
with the initial condition Q q (zi n i) = n 2 (l + Zi n i)~ 2 /4 at 
Zi n i = 2000 and substituting the resultant fl q (z) into (4), 
the corresponding E(z) can be obtained. 

For the NADE model, the cosmological constraints 
are found to be: n = 2.886^q o82 a ^ the lcr level and 
n = 2.886l°; 1 g3 at the 2cr lcvcL Thc bcst nt S ivcs val " 
ues of Xmm = 571.338, h = 0.685 and O to0 = 0.265. In 
Figure 1, we plot the relation of x 2 — n f° r the NADE 
model. 

According to [35], the NADE model is always consid- 
ered as a single-parameter model, and the initial con- 
dition is chosen to be Q, q (zi n i) = n 2 (l + Zi n i)~ 2 /A at 
Zi n i — 2000. We can, however, adopt the different 
perspective of thc NADE model being a two-parameter 
model. In doing so, we are interested in the result the 
observational data will tell us. To see this, we choose the 
initial condition il q o = 1 — f2 TO o for the NADE model, and 
then it becomes a two-parameter model with the free pa- 
rameters fi,„o and n. In Figure 2 (Left), we plot contours 
of the 68.3% and 95.4% confidence levels in the Q m o — n 
plane for the NADE model. The fit values for the model 
parameters are also shown in Table I. 

Now, with thc NADE model being viewed as a two- 
parameter model, it is of great interest to make a di- 
rect comparison with thc holographic dark energy (HDE) 
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FIG. 2: (color online) The probability contours at la and 2a confidence levels in the Q m o — n plane for the NADE model (Left) 
and in the Q m o — c plane for the HDE model (Right). Note that here the NADE model is regarded as a two-parameter model. 
The star denotes the best fit. 



TABLE I: The fit values for the NADE and HDE models. Note that here the NADE model is regarded as a two-parameter 
model. 



Model 








n/c 




Xmin 


NADE 


0.270™ (la 


4-U.U36 
-0.033 


(2a) 


206.762™" (la) 


-71.610 ( 2a > 


542.915 


HDE 


0.276± u ; U 02o (I" 


+0.036 

-0.033 


(2a) 


0.748± U U - (la)l 


-0.134 ( 2(J ) 


543.056 




FIG. 3: (color online) Left panel: the plot of x 2 versus n for the two-parameter NADE model with fi m0 = 0.270. Right panel: 
the plot of Q q versus z for the two-parameter NADE model with n varying from 137.702 to 337.974. 
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model. The equation of motion for the HDE fractional 
density Da is given by [30] 



n^ = n A (i-fi A ) 1 + -VU1 , (34) 



where c is a numerical parameter similar to n in the 
NADE model. Using ^ = - (1 + z) we get 

(35) 



— = -— — (1 - n A ) 1 + - v^a 

dz 1 + z \ c 



In Figure 2 (Right), we plot contours of the lcr and 2cr 
confidence levels in the fi m o — c plane for the HDE model. 
The fit values for the model parameters are also presented 
in Table I. 

Comparing with the HDE model, we find that the 
NADE model has a lower Xmin- Notwithstanding, it is 
obviously seen from Figure 3 (Left) that the data does not 
effectively constrain the parameter n, i.e., a very large 
range of values of n are allowed by the data. Figure 3 
(Left) shows the plot of x 2 versus n with fixed fl m o (which 
is fixed to the best fit value of 0.270). One should notice 
that x 2 tends to be a constant as n becomes large and 
the value of the constant \ 2 is just slightly larger than 
Xmin- ^ is, therefore, not surprising that n can range in 
value from 137.702 to 337.974 at the lcr level. Now, let 
us discuss the cosmological consequence of such a strong 
degeneracy of n. For this purpose, we plot the evolution 
of £l q (z) with n varying from 137.702 (the lower limit 
value at lcr) to 337.974 (the upper limit value at lcr) in 
Figure 3 (Right). We see that the curves with different 
n are almost totally degenerate in a narrow region. This 
indicates that the cosmological evolution tends to be the 
same when n takes large values in this model. In fact, 
we can also infer this from (7). We notice that the term 



2(1 + z) 



is negligible when n is large enough. Thus, 



(7) reduces to 



dCl 



1 



(36) 



Solving this equation we obtain p q = constant. There- 
fore, from the above analysis, we find that when the 
NADE model is regarded as a two-parameter model, the 
dark energy is more likely to behave as a cosmological 
constant. So, in the rest of this paper, we confine our 
discussions to the single-parameter NADE model. 

Next, we discuss the cases of the INADE model. Differ- 
ent cases of this model are denoted as INADE 1, INADE2, 
INADE3 and INADE4, respectively. Table II summa- 
rizes the fitting results for the four cases of the INADE 
model. For comparison, we also list the results of the 
NADE model. In this table we show the best fit, lcr and 
2cr values of the parameters and the Xmin values of the 
models. The best fit gives values of h = 0.692, 0.690, 
0.691 and 0.688 and n m0 = 0.240, 0.236, 0.237 and 0.239 
for the four interacting cases. One can see from Table II 
that the INADE 1, INADE2 and INADE3 models have a 



similar Xmin> an d tne INADE4 model gives a larger Xmin 
than the other three, but all are lower than the NADE 
model. In addition, a distinctive feature in the INADE4 
model is that the fit values of parameter b are all negative 
in the 2cr range. As discussed in the previous section, a 
negative b would lead to a negative p m in the future. So, 
in this sense, we have shown that the INADE4 model 
is not a reasonable model according to the observational 
data analysis. For the lcr and 2cr contours in the n — b 
plane for the four INADE models, we refer the reader 
to Figure 5. We will now discuss the high-z cosmic age 
problem caused by the old quasar APM 08279+5255 at 
redshift z = 3.91. 



IV. AGE PROBLEM: CHALLENGE AND WAY 
OUT 

The age of the universe at redshift z is given by 

dz' 



t(z) 



(37) 



(l + z')H(z') 

It is convenient to introduce a dimensionlcss cosmic age 



T cos (z) = H t(z) 



dz' 



(l + z')E(z')' 



(38) 



where E(z) = H(z)/H . At any given redshift, the age of 
the universe should not be less than the age of any object 
in the universe, namely, T cos (z) > T Q bj(z) = Hot bj(z), 
where t a bj(z) is the age of an OHRO at redshift z. For 
convenience, we define a dimensionlcss quantity, the ratio 
of the cosmic and OHRO ages, 



~ T ob 3 (z) 



')E(z') 



(39) 



Thus, the condition t(z) > 1 is equivalent to T cos (z) > 
Tobj(z). 

First, we will test the NADE model with the ages of 
the OHROs. We will use three OHROs: the old galaxy 
LBDS 53W091 at redshift z = 1.55, the old galaxy 
LBDS 53W069 at z = 1.43 and the old quasar APM 
08279+5255 at z = 3.91. The ages of the OHROs at 
z = 1.55 and z = 1.43 are 3.5 Gyr and 4.0 Gyr, respec- 
tively. For the age of the OHRO at z = 3.91, following 
[25], we use the lower bound estimate of t bj (3. 91) = 2.0 
Gyr. We calculate the age of the universe at differ- 
ent redshifts using the best fit results of the NADE 
model, n = 2.886 and h = 0.685, and then we obtain 
the values of r: r(1.55) = 1.215, r(1.43) = 1.136, and 
r(3.91) = 0.833, also shown in Table III. So, the NADE 
model can easily accommodate the OHROs at z = 1.55 
and 1.43, but cannot accommodate the old quasar at 
z = 3.91. Of course, the above result is only for the 
best fit. We will now see if the old quasar can be ac- 
commodated within the 2cr range. We show the result in 
Figure 4. In this figure, the blue line represents r(3.91) 
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TABLE II: The fit values for the NADE and INADE models. 



Model 




n 






b 






/Cmin 


NADE 


9 eoc+0.084 
^■»8t>_ .082 


(i^: 


,+U.lBU 
'-0.163 


(2a) 


N/A 






571.338 


INADE1 


3 1QQ+U.1M 




+0.324 

-0.290 


(2a) 


0.029™ (la) 


+0.014 

-0.015 


(2a) 


552.674 


INADE2 


q O/i £+0.218 
o.^io_„ 202 




.+0.364 

-0.322 


(2a) 


0.016™ (m: 


fU.UlU 
-0.010 


(2a) 


556.492 


INADE3 


q 9 q fi +0.207 
°- ZJU -0 .193 




,+0.344 

-0.307 


(2a) 


0.010™ (la) 


+0.006 

-0.006 


(2a) 


555.015 


INADE4 


3.208±«;^ 


(i^: 


+0.406 

'-0.339 


(2a) 


-0.027±^;^g (la 


,+0.022 fo \ 
'-0.025 ( 2<J ) 


561.634 



with n running over the 2a range; for simplicity, we fix 
h to the best fit value in the calculation. It is clear that 
the NADE model does not accommodate the old quasar 
APM 08279+5255, therefore, this age problem raises a 
serious challenge to the NADE model. 

To overcome this difficulty, we seek help from the pos- 
sible interaction between dark energy and dark matter. 
Thus, we next explore whether the old quasar APM 
08279+5255 can be accommodated in the INADE model. 
Based on the above results of observational constraints, 
we can easily accomplish this task. Using the best-fit 
values, we obtain the values of t(z) for the four INADE 
models, listed in Table III. One can see that the t(z) 
values indeed increase when the interaction is included 
in the model. Notwithstanding, for the old quasar at 
z = 3.91, the values of r(3.91) only increase from 0.83 
to about 0.92, still not exceeding the value 1. So, if 
only considering the best fit, the INADE models can- 
not solve the age crisis. Of course, it is evident that the 
age problem has been greatly alleviated via the interac- 
tion between dark energy and dark matter. As the next 
step, we will scan the full parameter space to explore 
whether there exists an area able to realize r(3.91) > 1. 
We show the result in Figure 5. The red line indicates 
where r(3.91) = 1, using the best fit value for h, and thus 
the region to the right denotes r(3.91) > 1, as indicated 
by the arrows. From this figure, one can clearly see that 
the red line passes through the 2a region of the n — b 
parameter space of the first three INADE models. For 
the fourth case, INADE4, the red line only intersects the 
edge of the 2a region. This result again shows that the 
form of interaction Q = 3bH(p q — p m ) is not reasonable, 
since it not only has a negative b, but also fails to provide 
a solution to the high-z cosmic age problem. 

The above analysis indicates that the old quasar APM 
08279+5255 can be successfully accommodated by the 
INADE model (at least for three cases) at the 2a level. 
There indeed exists an area within the 2a scope where 
r(3.91) > 1 is realized, and thus, it successfully solves 
the age crisis of the NADE model. To be modest, we do 
not assert that the cosmic age crisis has been completely 
solved by the INADE model. After all, only a small area 
and not the full region of the 2a scope resides to the right 
of the r(3.91) = 1 red line. However, it should be noted 
that the age problem has been significantly alleviated 
by the INADE model under the current observational 



constraints. 

Finally, we feel that it would be better if some ad- 
ditional comments on the age problem are given. The 
age problem is seen in the apparent old age of a single 
quasar, APM 08279+5255, whose high Fe/O ratio de- 
rived from X-ray observations has been used to argue for 
an age of at least 2 Gyr. At the quasar's redshift of 3.91, 
the standard ACDM model is only 1.6 Gyr old and thus 
apparently younger than the quasar itself. However, the 
age estimate of the quasar rests on the assumption that 
star formation and evolution proceeds in a way similar to 
that observed in the late universe. It is not even the ab- 
solute amount of Fe that serves as an argument, but the 
relative Fc/O abundance. Sub- millimeter observations 
show, however, that star formation seems to proceed in 
an atypical way in APM 08279+5255 [51]. This sheds 
doubt on the age estimate itself as well as on its relia- 
bility. Nevertheless, until a more accurate measurement 
of the age of the quasar is available, theorists have to 
seriously take the age problem into account. Likewise, 
a similar problem of star formation and evolution also 
exists for the high-z SNIa, but the observational result 
of the cosmic acceleration is still the mainstream. So, it 
is believed that the discussion of the age problem in the 
NADE model is meaningful and significant at the cur- 
rent stage. Of course, we do expect that future, more 
accurate measurements of the age of this old quasar will 
naturally eliminate the age crisis in dark energy models. 



V. CONCLUSION 

The agegraphic dark energy model originates from the 
holographic principle of quantum gravity, so it has phys- 
ical significance not only in a phenomenological aspect, 
but also in a theoretical one. Though the NADE model 
has been discussed extensively, it undoubtedly deserves 
further investigations. In this paper, we discuss the age 
problem in the NADE model. 

There is a lot of work addressing the age problem 
caused by the old quasar APM 08279+5255 at redshift 
z = 3.91, because this quasar has led many dark energy 
models into trouble, the ACDM and holographic dark 
energy models being no exception. We found in this pa- 
per that the NADE model is also afflicted by the age 
problem. So, we explore whether the inclusion of the in- 
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FIG. 4: (color online) The plot of t(3.91) versus n for the NADE model. The region avoiding the age problem is that above 
the t(3.91) = 1 line (red). 



TABLE III: The values of t(z) in the NADE, INADE1, INADE2, INADE3 and INADE4 models at best fit concerning LBDS 
53W091 at z = 1.55, LBDS 53W069 at z = 1.43, and the old quasar APM 08279+5255 at z = 3.91. 



Model (n, b, h) 


r(1.55) 


r (1.43) 


r (3.91) 


NADE (2.886, 0, 0.685) 


1.215 


1.136 


0.833 


INADE1 (3.199, 0.029, 0.692) 


1.326 


1.238 


0.919 


INADE2 (3.245, 0.016, 0.690) 


1.318 


1.230 


0.923 


INADE3 (3.236, 0.010, 0.691) 


1.323 


1.235 


0.923 



INADE4 (3.208, -0.027, 0.688) 1.290 1.204 0.908 



teraction between dark energy and dark matter in this 
model can solve the age problem. 

First, we used current observational data to constrain 
the NADE and INADE models. For the data, we use the 
SNIa Union2 sample, the CMB shift parameter R from 
7-yr WMAP, and the BAO parameter A from the SDSS. 
We determined the age of the universe in the NADE 
model by using the fit results and found that the NADE 
model cannot realize r(3.91) > 1 within the full parame- 
ter space. We also explored the cosmological consequence 
of regarding the NADE model as a two-parameter model. 
The fit results show that in such a model the dark energy 
is more likely to behave as a cosmological constant. 

For the possible interaction between dark energy and 
dark matter, we considered four phcnomenological cases: 
Q = 3bHp q , Q = 3bHp m , Q = 3bH(p q + p m ) and 
Q = 3bH{p q — p m ). The observational data favor a pos- 
itive b for the three cases, but give a negative b for the 
fourth case. We found that when the interaction is taken 
into account, the age problem caused by the quasar can 
be successfully solved in the first three cases of the IN- 
ADE model, at the 2a level. The isoline r(3.91) = 1 
passes through the 2a region of the n — b parameter 



space of the INADE model. So, there indeed exists an 
area within the 2a scope realizing t(3.91) > 1 and thus 
this model successfully solves the age crisis of the NADE 
model. Though we cannot assert that the age crisis has 
been eliminated in the INADE model, since most of the 
parameter space is still in the troubled area and only 
the lower bound age of the quasar is used in the test, 
we are sure that the age problem has been significantly 
alleviated in the INADE model under the current obser- 
vational constraints. Therefore, our work can be viewed 
as further support of the INADE model. 
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